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The thermal fluctations of vortices in a superconductor can be usefuUy mapped onto the quantum 
fluctations of a collection of bosons at T = moving in 2 dimensions. When the superconductor is 
a thin platelet with the magnetic field parallel to its surface, the interacting quantum bosons are 
effectively moving in 1 dimension, allowing for powerful Luttinger liquid methods to be applied. 
Here we consider how this 1 dimensional limit is approached, studying the interaction of vortices 
with the platelet surfaces and each other. Using realistic parameters and vortex interactions for an 
' underdoped YBCO platelet we determine the scattering length, a, characterizing the low energy 

interaction of a vortex pair as a function of the platelet thickness, a determines the Luttinger 
, parameter, g, for the quantum system at low densities, no: g — > 1 — 2ano. 
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PACS numbers: 72.15.Qm, 73.21. La, 73.23.Hk 



^ I. INTRODUCTION 

o, 

O , Thermal fluctuations of vortices, taking into account pinning by impurities and vortex-vortex interactions, is a 
^ • challenging and technologically important problem in statistical physics. An elegant approach to this subject is to 
Oh' map each fluctating vortex line into the world line of a quantum particle in a Feynman path integral, with the magnetic 
^ . field direction becoming the imaginary time direction and the bosons moving in the other two spatial directionsJ^ 
^ ' Such a mapping is especially powerful for studying columnar defects which become static point defects in the quantum 
"j^ , model. When the field direction is tilted relative to the (parallel) pins a novel non-Hermitean 2-dimensional many- 
' body quantum problem arises If the superconductor is a thin platelet, of thickness of order the penetration depth, 
I , with the field lying in the plane, then the quantum bosons are essentially restricted to one dimension. This allows 
' theoretical techniques including Tomonoga Luttinger liquid (TLL) theory and Density Matrix Renormalization Group 
Ch [ (DMRG) to be brought to bear^ rendering tractible a formidable problem. It may be feasible to realize this classical 
analogue of a Luttinger liquid experimentally using high-Tc superconductors. A promising candidate would be a very 
clean highly underdoped YBCO single crystal in which the penetration depth Ac can be as large as 50 microns. A 
platelet should be cleaved with thickness in the a-direction of order .1 to 1 mm. and then a magnetic field should be 
applied in the fe-direction. 

^ , It was shown in ^jSj] that critical phenomena connected with rotating the field direction away from the pin direction 
is controlled by the dimensionlcss Luttinger parameter, g. In the case g > 1 columnar defects are irrelevant and have 
little effect on the long-distance properties of the vortices. On the other hand, for g < 1 they are relevant and an 
i arbitrarily weak pinning potential drastically alters the system. Thus it is of considerable interest to determine g and 
how it depends on the parameters of the system, including the density, hq. In the dilute limit g approaches unity 
and the interacting boson system becomes equivalent to non-interacting fcrmions. The leading density dependent 



g = l-2ano. (1.1) 

Here a is the 1-dimensional scattering length. This is defined in terms of the 1-dimcnsional even-channel phase shift, 
S{k). The even wave functions have the asymptotic long-distance behaviour: 

^eix) ^ sm[k\x\ ~ 6{k)], (1.2) 

where x is the separation of the 2 bosons. At fc ^ the phase shift is linear in k: 

5{k) ak, (1.3) 

implying: 

^pe{x) smk{\x\ — a). (1.4) 

Thus the relevance or irrelevance of pinning, in the dilute limit is determined by the sign of a. It is important to realize 
that this crucial sign is not fixed by the requirement that the boson-boson interaction be repulsive. For example, 
an infinite hard core repulsion of range aq leads to a scattering length a = oq > 0. On the other hand, a repulsive 
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(5- function interaction, v5{x), leads to a negative scattering length, a = —l/{fiv) (where fj, is the reduced mass). In a 
confined geometry the 1 dimensional scattering length depends not only on the direct inter-particle interaction but 
also on the effects of the boundaries. This problem was solved by Olshanii^ for the case of ultra-cold atoms in a 
harmonic cylindrical trap, where it was shown that the sign of a can be positive or negative depending on the ratio 
of the (positive) 3 dimensional scattering length to the trap radius. 

In this paper we study the properties of two interacting vortices in a thin platelet, or equivalently of 2 interacting 
bosons restricted to a narrow strip. We begin with the usual modified Bessel-function interaction between vortices 
given by anisotropic London theory. Standard boundary conditions at the edges of the platelet imply the existence 
of an infinite set of image vortices for each physical vortex. The interaction of an isolated vortex with its images, and 
with the external magnetic field determines its wave- function, /(y), in the quantum mechanical analogue, determining 
the probability of the vortex being at a distance y from the centre of the platelet. It is peaked near the centre, y = 0- 
We then consider the scattering of two physical vortices, taking into account the interactions with all image vortices. 
In general the two vortices could move off centre (away from yi = y2 = 0) as they scatter. Thus the calculation of 
the effective ID scattering length requires, in principle, solving for a 2-dimensional 2-body wave-function. 

Fortunately, there is a very large dimensionless number that appears quite generally in the thermodynamics of 
vortices, and which simplifies our calculations considerably. Consider, for simplicity, a macroscopic isotropic London 
superconductor of penetration depth A and coherence length ^. We approximate the Gibbs free energy for N vortices 
as: 



(1.5) 



Here (po = hc/{2e) « 2 x 10 ^G.-cm.^ is the flux quantum and 



the tilt modulus, is simply the energy per unit length of the vortexii r is the spatial co-ordinate along the field 
direction and ri{T) describes the shape of the i*'' vortex. We have approximated the vortex- vortex interaction as only 
depending on the difference of the r^'s at the same value of r and used the standard London model result for the 
interaction energy per unit length between straight parallel vorticesj^ given by the modified Bessel function, Kq. (This 
needs to be cut off at short distances of order w ^.) We approximate the partition function by an integral over 
vortex paths, ri(r), weighted by the Boltzmann factor, exp[— G/(A:bT)]. By identifying G/{kBT) with S/h where 5" is 
the classical action for N interacting bosons, the classical partition function describing thermal fluctuations of vortices 
becomes equivalent to the Feynman path integral for interacting bosons. In this way various thermal properties of the 
vortex system can be conveniently obtained from the quantum system»i— (We set h and ks = 1.) The corresponding 
Hamiltonian is: 

1 ^ dy^ 



2m ^ ' Stt^A^T 

2 — 1 2<_7 



with 



h/T. (1.8) 



(Strictly speaking, even after inserting appropriate factors of h and fc^, both terms in this Hamiltonian have dimensions 
of inverse length rather than energy. This can be traced back to the fact that in Eq. (|1.5p r is a spatial coordinate in 
the classical model but is treated as an imaginary time in the analogue quantum one. This creates no problems for 
our analysis since physical quantities calculated using this quantum approach involve appropriate ratios of parameters 
with the correction dimensions, as we shall see.) It is convenient to change to dimensionless length variables, letting: 

u., = (1.9) 

We then may write the Hamiltonian in dimensionless form: 

N 

2mX^H = -^\/^^+Vo^Ko{\u,~Uj\) (1.10) 

2—1 2<j 
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FIG. 1: A thin platelet superconductor. 



where the dimensionless parameter which measures the interaction strength is: 



Noting that 

(t>l/{?>-K^kB) = 3.9223 X lO'^ ^wl - K (1.12) 

we see that Vq 3> 1 for essentially any superconductor at any T < T^. This means that the analogue quantum 
mechanical bosons have very strong short-range interactions when measured in dimensionless units. Variants of this 
large number will appear when we consider the potential energy function that holds the vortices in the middle of the 
slab and the interaction between vortices inside the slab. This implies that the vortices stay near the centre of the 
slab up to rather large slab widths justifying a 1 dimensional approximaton. It also allows an unusual but powerful 
semi-classical approximation to be applied to the 1 dimensional problem yielding an explicit formula for the scattering 
length as a function the platelet thickness d and other parameters (Aa, Ac, ^a, and T). Our conclusion is that a is 
positive and large for narrow platelets, increasing with d and having a value a w 19Aa for d ~ lOXc ~ .5mm. This 
implies that columnar pins arc relevant and also that the system rapidly leaves the dilute regime at low densities of 
order l/(20Aa), corresponding to fields, H only slightly above Hd. 

In the next section we briefly review London theory and discuss properties of a single vortex in a thin platelet. In 
Sec. Ill we consider 2 interacting vortices, determining the scattering length. Sec. IV contains conclusions. 



II. A SINGLE VORTEX IN A THIN PLATELET SUPERCONDUCTOR 



We make the London approximatioufi valid when the penetration depth is much longer than the coherence length, 
A ^ ^. We label the direction perpendicular to the platelet the y-direction, and label the direction of the magnetic 
field the z direction. In a YBCO crystal the most promising geometry may be choosing y and z to be the a and b 
directions (or vice versa). See Fig. ([1]). Thus the thin direction of the platelet is the a direction, not the usual growth 
direction, which is c. Such a sample could presumably be obtained by cleaving a macroscopic sample. The magnetic 
field of a single vortex, centered at r = {x, y) = thus obeys: 

^^oS^r). (2.1) 

The Dirac (5-function at the vortex core should actually be smeared over a distance of order ^, the coherence length. 
Note that the decay of the magnetic field in the x = c direction is governed by supercurrcnts running in the in the 
y = a direction and hence involves A^ whereas the decay in the y = a direction is governed by supercurrcnts running 
in the x ~ c direction and hence involves Ac- In extremely underdopcd YBCO crystals typical parameter values arc 

Ac = 50/-im 

Xa = .5//m 

= 5nm 

^c = .05nm 

Tc = 17K. (2.2) 
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(The value of may be a bit small compared to existing measurements but it is convenient to assume the result 
which follows from anisotropic Ginsburg-Landau theory: Ac/Aq = Ca/Cc- In any event our results only depend 
logarithmically on the ^'s.) Thus the vortex is extremely elliptical: much more extended in the y = a direction. To 
acheive the two-dimensional limit, we need the sample thickness to be of order the vortex size. (Actually, as we shall 
see a thickness of up to ten times the vortex size or more is alright.) Thus we can take advantage of the larger Ac 
by cleaving our crystal in the a-direction. We will refer to the parameters in Eq. (j2.2p at temperature T fa Tc, with 
a = y the thin direction as the standard parameters. However, our results should also apply at lower temperatures. 
Note that the vortex wandering that we arc concerned with here occurs primarily in the x = c direction. The vortices 
presumably feel a periodic potential with wave-length given by the lattice constant in the c-dircction. We will ignore 
this here. In the dilute limit that we are considering it is not expected to have an important effect. 

The solution of Eq. (|2.ip at distances r ^ ^, ignoring for now the boundaries, can be found by using a Fourier 
transform and is 



47r2 



OO pOO 



oo ^ — cxD '*a x 




(2.3) 
(2.4) 



where Kq is a modified Bcsscl function. 
The energy per unit length of a vortex is 



da. 



where the integral is taken over an ellipse of radii ^a, £,c around the vortex core. Thus 

'A. 



167r2AaA. 



■In 



(2.5) 



(2.6) 



to logarithmic accuracy. For the parameters in Eq. (|2.2[) we have e = 10~^erg/cm. 

We remark that the interaction energy per unit length between 2 straight parallel vortices separated by a vector r 
is simply: 



47r 



(2.7) 



where hi2{'r) is the magnetic field at the location of one vortex produced by the other, Eq. (|2.4p . 

We now consider a single straight vortex in an infinite slab, of thickness d, extending from —d/2 to d/2. The 
presence of boundaries of the superconductor at y = ±d/2, imposes boundary conditions. 



= Jy 



dh 
dx 



(2.8) 



A vortex at position j/ in a slab of thickness d creates image vortices with magnetic field in the opposite direction at 
positions (2n -I- l)d — y for integral n, and creates image vortices with field in the same direction at positions 2nd + y 
for all nonzero integral n. Thus the magnetic field at location (a;, 7/2) for a vortex at (0,j/i) is: 



h{x;y2,yi) 



27rAaAc 



71 — — OQ 



x\'^ f 2nd + j/i - 2/2 ^ ^ 



x \ ^ ( (2n -I- l)c? — j/i — 2/2 ^ 



A. 



(2.9) 



In addition to the field of the vortex and its images an additional field occurs inside the superconductor when a field 
H is applied outside of it: 



(2.10) 
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The magnetic field at the position of the vortex, (0, y) due to all its image vortices is: 



h2{y) 



27rAaAc 



E ^0 



ri=±l,±3... 

The Gibbs free energy depends on the position y of the vortex as: 



n=±2,±4... 



A. 



^i(y) = -r-{hi{y) + ^^2(2/)] + constant. 
47r z 



(2.11) 



(2.12) 



This is plotted in Fig. (2(a)) and (2(b) at H = Hd- It has a minimum at y — 0, large barriers at intermediate y 
and then appears to diverge to — oo at y ^ ±d/2: 



167r2AaA, 



■ln[Ac/(d±2y)]. 



(2.13) 



This divergence is due to the interaction of the vortex with its image at y^d. This divergence should actually be cut 
off at y^d/2 of order ^ due to corrections to London theory. We take this into account by replacing K(j[{nd — 2y) / \^ 
by KQ[{nd — 2y)/X^ — KQ[{nd — 2y)/^] in Eq. (|2.1ip . For H > H^i, we might expect the true minimum energy to be 
at y = 0. This gives the usual formula for Hci for an anistropic superconductor: 



tlcl ~ — £• 
00 



(2.14) 



where e is the energy per unit length of a bulk vortex, Eq. (|2.6p . 
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FIG. 2: (a) and (b):The Gibbs potential, Eqs. (|2.10|) - (|2.12|l . for a single vortex with our standard parameters, H = Hci and 
d = 10, lOOAc; y is in units of Ac. (c) and (d): The corresponding wave functions ,f{y). 

We now turn to a study of thermal fluctuations of a single vortex inside the slab for such an anisotropic supercon- 
ductor. With the magnetic field along the 6-axis, the tilt modulus is very different for vortex tilting in the a = y or 
c = X direction. The elastic energy is written: 



Go 



dT 



2 



(2.15) 



Due to the assumed symmetry under rotations in the a-b plane, the tilt modulus for tilting in the a direction is just 
given by the energy per unit length: 



(2.16) 
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where e is given in Eq. (|2.6p . On the other hand, the energy per unit length of a vortex ahgned parallel to the c = x 
axis is much larger, resulting in the tilt modulus^ 



(2.17) 



To this must be added the y-dependent free energy: 



Gi = / drVMr)]- 



(2.18) 

Here we have again considered only "instantaneous" interactions between the vortex and its images, at a fixed value 
of T. We do not expect this approximation to qualitatively change the long distance physics in the dilute limit. If we 
consider a very long vortex, in a sample of macroscopic length in the z = t = b direction, then the probability of the 
displacement of the vortex from the centre of the slab having some value y is simply given by |/(y)P where / is the 
ground state wave-function of the one-dimensional Hamiltonian: 



Hi=- 



1 



2ma 



Vijy) 

T 



with 



ma = —. 



It is now convenient to define the dimensionless length variable: 



in terms of which: 



y = y/K, 



^o[(2n + l)d/Ac - 2y] + HK/ia) 



cosh(y/2) 
cosh((i/2Ac~ 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



where 



167r2Aar2 



IGTT^AaT 



In 



(2.23) 



[We have set H — Hd given in Eq. (|2.14p .] Since we will choose d of order Ac, we see that the dimensionless number 
Vbj, characterizes the height of the barriers holding the vortex at the center of the platelet. Using the numbers in Eq. 
and choosing T = Tc we find: 



Vr 



Oy 



(2.24) 



We solve this Schroedingcr equation numerically for the groundstatc, for d = lOAc, Fig. [2(c)| , finding that the particle 
makes only very small quantum fluctuations away from y = due to the huge barriers. For the above parameters we 
find 



> 



.01419Ap = .001419d. 



(2.25) 



We should estimate the true critical field, Hci{d) using the ground state energy of the quantum Hamiltonian. If the 
field is too low the particle can tunnel through the barrier corresponding to a vortex terminating at some value of 
r with the magnetic flux leaving the superconductor. However, due to the exceedingly high barrier, the tunnelling 
probability is miniscule and the system will not be very sensitive to the precise value of H . We should also remark 
that, to solve the Schroedinger equation precisely we need to specifiy some boundary conditions on the wave- function 
at y = ±d/2. We imposed vanishing boundary conditions. Fortunately, the extremely large dimensionless barrier also 
renders our results very insensitive to this choice. 
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III. TWO VORTICES 

Consider two straight parallel vortices inside the platelet, at locations {xi,yi). The Gibbs free energy per unit 
length is: 

V{xi - X2;yi,y2) = -^^^i^i - X2;yi,y2) + Vi(yi) + 1^2(2/2)- (3.1) 

By translational invariance in the x-direction, V depends only on the difference of the x-coordinates of the two 
vortices: 

X = xi — X2- (3.2) 

Here h{x,yi,y2) is given in Eq. (|2.9|) and Vi{y) by Eq. (|2.12p . The first term in Eq. ()3.ip represents the interaction 
of one vortex with the other one and with the images of the other one. The second and third terms represent the 
interaction of each vortex with its own images and with the screened external field. V{x, yi, 2/2) has a deep minimum 
at ?/i = ?/2 = for large |x| ^ Aa- V^(x-,0,0) has a large peak centred at a; = with a weak, logarithmic divergence 
right at a; = 0. This logarithmic divergence should be cut off at scales of order ^; however this cut off has essentially 
no effect on the scattering length, as we shall see. 

Again we may study the thermodynamics of two wiggling vortices in the platelet by mapping onto a quantum 
mechanics model. We make the fundamental assumption that the vortices only bend on long length scales (compared 
to Aa) and that we may approximate the vortex- vortex (and vortex-image vortex) interaction by an "instantaneous" 
one at a fixed value of r = z. The Boltzmann sum is now over the configuration of two vortices and we must include 
the vortex-vortex interaction in the free energy. Again identifying the free energy with the imaginary time action, we 
see that the corresponding quantum Hamiltonian is: 



H = - 

Here rua is given in Eq. (|2.20[) and 



1 ^ 

4=1 



1 / y I f d 

ma V dyi J mc V '^^i 



nx^^_ (3.3) 



"^c^l (3.4) 



where is given in Eq. (|2.17p . In the quantum analogue, the vortices obey Bose statistics^ and consequently 
the two-body wave-function must be symmetric: even under x —x. The asymptotic behavior of the low energy 
wave-functions at 3> Aa is given by: 

i^{x, yi,y2) /(2/i)/(2/2) sin[fc|2;| - S{k)] (3.5) 

where /(y) is the ground state wave-function for a single vortex, discussed in Sec. H. The energy of this scattering 
state is: 

E = 2E, + ^^ (3.6) 

where Ei is the ground state energy for a single vortex, discussed in the previous section and the reduced mass which 
governs the relative motion is 

= m,/2. (3.7) 

The scattering length is defined by Eq. (|1.3p . 

It turns out that, due to the large barrier near a: = 0, for small yi, the scattering length is determined almost 
completely by the large x asymptotic form of the potential, until d gets very large compared to Ac- The x-dependcnt 
part of the exact potential can be written as a sum of exponentials: 

V{x; yi, 2/2) = f T' f ) , , ;,^2/^ + ^^^^^^^ + ^1(2/2). (3.8) 

where we define 

fin, y,/d, y2/d) ^ { cos{nnyJd)cos(^nny2/d) if n is odd 

J \ ■, ui-i ■, a-^i I 1 sm(n7ryi/d) sm(n7ry2/") if n is even ^ ' 
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At large |a;| we may approximate the sum by the first term only: 



where we have defined, for convenience, a reduced value of Xa- 



-1-1/2 



(3.10) 



(3.11) 



For this approximation to hold, we need that the first term dominates all other terms. In the one-dimensional limit 
(i.e. 2/1=2/2 = 0) the condition for large x is 



exp 



> 1. 



(3.12) 



Note that, unlike the direct vortex-vortex interaction, this potential has a simple exponential dependence on x at 
large |x| albeit with a reduced penetration depth. 

Due to the large barriers in Vi{y) we expect the low energy scattering states to be confined to j/i w 0. We first 
calculate a assuming that the vortices stay exactly in the middle of the slab, iji = throughout the scattering process. 
We return to a further discussion of why this is reasonable at the end of this section. This reduces the problem to a 
one-dimcnsional quantum mechanics model with Hamiltonian: 



H 



1 (P V{x) 
' 2n dx^ ^ T 



where 



V{x) = V{x;0,0) 



^-\x\/X 



(3.13) 



(3.14) 



We look for parity even solutions of this Schroedinger equation with asymptotic behavior il}{x) — > sin[A;|a;| — 5{k)\ 
with 5{k) ak ds k ^ Q. Note that in the small k limit, ^'(2^) ^ sin[fc(|x| — a)] for x ':$> X. Then, if we consider an 
intermediate range of x: 



A < < 1/fc, 
we may approximate the wave-function by a linear form: 

ipix) (X \x\ — a. 

Thus, to find the scattering length we need to simply solve the zero energy Schroedinger equation: 

1 d2 ^(3.)- 



T 



t/j = 0. 



(3.15) 
(3.16) 

(3.17) 



Note that this reduces the eigenvalue problem to a simple initial value problem. We simply impose the initial 
conditions: 



dip 



(0) 







dx 

m = 1 



(3.18) 



and solve the zero energy equation. The asymptotic behavior of the solution at |a;| 3> A is given by Eq. p.l6p which 
determines the scattering length, a. 

It is convenient to introduce a dimcnsionlcss length variable: 



X = x/X 

in terms of which the Schroedinger equation becomes, at large |i|: 



i:{x) = 2^iX^EiP{i). 



(3.19) 



(3.20) 
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Here the dimcnsionlcss number which measures the strength of the repulsive potential is: 



Vr 



Ox 



1 / 



Xc/d 



[l + (7rA,/d)2]3/2- 



AnXldT^ 27r \ATr XaT ^ 
Using our estimates of the parameters in Eq. (|2.2|1 with T = Tc and d = lOAc we find: 

Vo^ = 1.07 X 10^ 



(3.21) 



(3.22) 



Note that Vqx oc 1/d aX d ^ Ac. Importantly Vqx S> 1 when d is of order Ac and remains large out to extremely 
large values of d/Xc- The largeness of Vqx leads to a large scattering length, allows for an unusual semi-classical 
solution approximation and also helps to justify setting y,; = as we shall see below. It is of course, possible to solve 
the Schroedinger equation numerically for specified values of the parameters. However, the largeness of V^x and Vby 
creates numerical difficulties for standard algorithms, when one attempts to solve the full 3-dimensional problem, 
including the yi. In this case it is much easier, and more transparent, to use the semi-classical approximation. 

Our semi-classical approximation for the scattering length at Vqx ^ 1 begins with the observation that the classical 
turning point for 2 particles approaching each other with a small relative momentum, k occurs at i 3> 1. Therefore 
a is determined almost completely by the large x form of the potential in Eq. (|3.10p . Using the large x form of the 
potential, we may solve the one-dimensional Schroedinger equation exactly. To do this we change variables to: 



u = 2v/^e-l*l/2. 
The zero energy Schroedinger equation, (j3.17[) . simplifies to 

u^ijj" + uijj — u^ijj = 



(3.23) 
(3.24) 



where the primes denote differentiation with respect to u. This is the zeroth order modified Bessel's differential 
equation, and the general solution is given in terms of the modified Bessel functions: 



■)p = ci/o(u) + C2Ko{u) 

= ci/o(2/^;:e-l^l/2^) + C2Xo(2yi^e-l^l/2^). 



For large u (i.e. x << InVp) we have 



'2'KU 



C27re 



or writing in terms of the original variables and putting back factors of A, we have 



1 



cie 



C27re 



On the other hand, for small u (i.e. x » InVo) we have 

« ci - C2(ln(u/2) + 7) = ci -f C2 - i In Vqx - 7 

where 7 w 0.5772 is the Euler-Mascheroni constant. Note that the last formula can be written as 

V'- ||(|a;| - a) 
where the scattering length a is (putting back factors of A) 

a = A(lnUo:r + 27-2ci/c2). 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



To determine Ci/c2 we match our solution in the region 1 <C |5:| <C InVox to the WKB solution which works for |a;| 
not too large, in the region where V{x) is large. The even WKB wave function for E' = is given by 



ipix) = A < exp 



^yV'{x')dx' 



+ exp 


[-1 







y/V'{x')dx' 



(3.31) 
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where 



V'{x) = 2fiV{x;0,0) 



(3.32) 



and V is the exact potential of Eq. (|3.ip . (Note that we do not make any large x approximation to V now.) We can 
rewrite this as 



ipix) = A < e " cxp 



y'V'{x')dx' 



e exp 



y/V'{x')dx' 



where 



^/V'{x)dx. 



(3.33) 



(3.34) 



Now, if X is large enough for our asymptotic expression V'{x) « VoxG j}? to hold, then the integral can be done 
quite readily: 



xl2\ 



ip{x) — A^e " exp 
Comparison with Eq. (|3.27|) thus gives 

Ci/c2 ~ 7re~^" — Trexp 



e exp 



'x/2\ 



y/V'{x)dx 



(3.35) 



(3.36) 



This quantity is exponentially small in the large quantity Vqx so it is completely negligible. Note also that the 
logarithmic divergence of V{x;Q,0) at a; has no important effects, leaving the integral finite in Eq. p.36p . The 
last term in Eq. (j3.30p essentiallv vanishes, and we simply have that 



a = XiliiVox + 2-f). 



(3.37) 



Asymptotically, the scattering length is linearly dependent on the logarithm of the size of the potential. 

Interestingly, we get almost the same result for the odd wave functions, except that the sign of Ci/c2 is reversed. The 
even channel and odd channel scattering lengths are therefore almost exactly the same. However, it is the difference 
between the even and odd channel scattering lengths that determines the transmission coefficient, and it is only then 
that Ci/c2 plays an important role. 

We have based this approximation on the assumption that there exists a region of separation x, such that the 
approximation Eq. (|3.27[) holds, and the WKB approximation to the wave function also holds. Typically the matching 
is done around the region a; « a; therefore we need that (using Eq. (|3.37p and Eq. p.l2p . and noticing A « Aq for 
large d) 



exp 



27rA, 



ln(e2^yo.) 



> 1 



(3.38) 



This is the condition that must be satisfied for the formula Eq. (|3.37p to hold. 

The wave-function calculated numerically to high precision (in the 1-dimensional approximation) and the semi- 
classical wave function are compared in Fig. ([3]) for our standard parameters and d = IOAq.As can be seen, the two 
wave functions give good agreement in the large- a; regime (with an error < 1% for a; > 18Aa). The semi-classical 
wave function is grossly inaccurate in the small- a; region (x < 15Aa) but the wave-function is neglegible there anyway. 
The predicted semi-classical scattering length is a = A(lnVb -I- 2j) = 18.7398Aa. The actual scattering length of 
the numerically determined wave function (obtained by fitting the wave function in the large- a; regime to a linear 
function) is 18.7409Aa. The semi-classical wave function gives an error of less than 0.01%. In Fig. ^ we show the 
scattering length versus d, comparing our numerical results to the semi-classical approximation (in both cases making 
the 1-dimensional approximation). Most of the c?-dependence in our semi-classical formula, Eq. (|3.37p . arises from 
the d-dependence of the reduced penetration depth, A, given in Eq. (|3.1ip . As the sample thickness decreases, the 
effective range of the interaction potential. A, also decreases, and the scattering length simply scales with it, up to 
logarithmic corrections coming from Vqx, defined in Eq. p.2ip . The semi-classical and numerical values agree within 
1% up to about d = 22Ac. For d = 22Ac, we have that (refer to Condition p.38p ) 



27rA, 



ln(e2''yo..) = 1.56 



(3.39) 
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FIG. 3: The wave functions in the k —> limit comparing a precise numerical calculation (in the 1 dimensional approximation) 
to the semi-classical approximation for our standard parameters and d — lOAc. Lengths are in units of Aa. 
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FIG. 4: The scattering lengths based on a precise numerical calculation (in the 1 dimensional approximation) compared to our 
semi-classical approximation for our standard parameters, d is in units of Ac, the scattering length in units of A^. 



which is already not far from unity. There will therefore be significant deviations of the true value from Eq. (j3.37p . We 
could also notice that the scattering length tends towards a finite value for large d. This is because as the thickness of 
the sample grows, the image vortices move farther away from the orginal vortices until their effects become negligible. 

Next, we discuss the validity of our 1-dimensional approximation, setting yi = 0, which is justified by the fact that 
the single vortex wave-function, f{y) is so sharply peaked near y — 0. If we look at the shape of the potentials, 
Vi{y), we see that for small d (Fig. (2(a) )) the potential has an obvious minimum in the center, and is approximately 
simple harmonic near the center. For larger d (Fig. (2(b))), however, the potential is almost negligible except for 
a large potential barrier close to (but not at) the edges; the potential will be qualitatively more similar to a square 
well. Therefore as d increases, we would expect the shape of the wave function to morph from a confined Gaussian 
to a spread-out sinusoidal (Fig. (2(c)). Fig. (2(d))). The spread of the single-vortex wave function {{y"^)^^^) is 
plotted in Fig. ([5]). The thickness reaches 1% of the platelet thickness at around d = 35Ac, an indication that our 
one-dimensional approximation fails above this value. The fact that the thickness becomes linearly dependent on d 
at large d also suggests that the wave function tends to a fixed shape (a sinusoidal). 

A more systematic approximation to solving the full 2-body 2-dimensional Schroedinger equation of Eq. 
would be to write: 



■<P{x;yi,y2) « tpiix) f (yi) f {y2) 



(3.40) 



where ipii^) is the 1-dimensional wave- function found above and /(j/i) is the single vortex wave- function. We could 
then improve our estimate of the 1-dimcnsional effective potential by using: 



V{x)^ / dyidy2\fiyi)Wiy2)\'V{x;yi,y2) 



(3.41) 



rather than simply V{x) w V{x;Q,0). However, because f{y) is so sharply peaked at ?/ w this makes a negligible 
difference. 
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FIG. 5: The spread of the single-vortex wave function, ^/<j/^>, plotted against d using standard parameters. Both lengths 
are in units of Ac. 



Note that our calculation of a depended essentially only on V{x) in the large x region, x ^ X. Our consideration of 
the small x region only served to determined ci/c2 which was exponentially small anyway and can simply be ignored. 
In this large x region, Eq. (|3.10p is a good approximation, the wave-function approximately factorizes and the large 
barriers in Vi{y) ensure that the wave-function is strongly peaked near j/; ~ 0. For smaller values of x the wave- 
function presumably spreads out more in the y direction. However, at smaller x the wave-function is exponentially 
small anyway. 

Finally, we consider the case where the thin direction of the YBCO platelet is the c direction: y=c (and x=a, z=b). 
In this case the roles of Xa and Ac are switched, as are the roles of and ic and we get 



Wtt^T^X^ 



1 



AirXldT^ 47r V 47rAaT 



Wn^XaT 



Xa/d 



[1 + (^AJci)2]3/2 



(3.42) 



where we now define: 



A = Ac 



ttAq 

d 



-1-1/2 



(3.43) 



Apart from some unimportant factors of 2, the formulas for Voy and Vqx are the same as for the other geometry except 
that d now appears in the dimensionless ratio Xa/d rather than Xc/d. So we now conclude that the 1 dimensional 
approximation holds for d < 20Xa ~ lOyum. and the semi-classical approximation holds out to roughly the same value 
of d. 



IV. CONCLUSIONS 



Our main result is the formula Eq. p.37p for the scattering length. This is plotted versus the platelet thickness, 
d, for our standard parmeters, in Fig. ([4]). Note that a is everywhere positive and a/Xa is everywhere quite large, 
having the value a/Xa ~ 19 at c? w lOAc. At somewhat larger values of d we expect our semi-classical approximation 
to the one-dimensional problem to break down and, more problematically, the one-dimensional approximation itself 
to start to fail. 

The main use of our formula for a is not, of course, for studying the system with only 2 vortices, but rather for 
studying the thermodynamic limit of many vortices. In the dilute limit, nga <C 1 (where uq is the vortex density 
per unit length), a determines the Luttinger parameter via Eq. p.ip . Of course the Luttinger liquid treatment of 
the problem assumes that it is fundamentally 1-dimensional. Our calculations here indicate that the 1-dimensional 
approximation should be good, at least in the dilute limit nga ^ 1, up to platelet thicknesses of order d = 10 Ac or 
more, since the vortices stay very close to the centre of the platelet. Furthermore, we have determined the Luttinger 
parameter for this range of thicknesses and vortex densities. When a is large the Luttinger parameter decreases 
rapidly for increasing vortex density. It was argued in HI) that, at high densites, g <^ 1. Taken together, these results 
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suggest a rapid monotonic drop of g from 1 with increasing density. In this case, columnar pins would be highly 
relevant for essentially all fields above iJ^i • Thus a promising region to look at experimentally might be very close to 
Hci with low vortex densities, uq ^ 1/Aa and samples of thickness around lOAc- 
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